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METHOD OF CAIiCULAT I NG THE FAST FOURIER TRANSFORM AND 
THE INVERSE FAST FOURIER TRANSFORM 



This invention relates to a method of calculating 
the fast Fourier transform or the inverse fast Fourier 
transform of a series of real numbers or a series of 
5 > conjugated complex samples . 

The Fourier transform is probably one of the most 
important tools for analyzing, designing and 
implementing signal processing algorithms, and the 
existence of efficient algorithms, such as that of the 
0 fast Fourier transform, has been a major factor for 
this situation. Although most Fourier transform 
algorithms are designed for transforming series of 
complex numbers, there are, however, various 
applications, such as image or acoustic signal 
5 processing or certain types of multicarrier modulation 
wherein the series to be transformed are real numbers. 

In general, the direct Fourier transform and the 
inverse Fourier transform respectively set up the 
following relations between two series of N complex 
0 numbers, x(n) and X (n) : 

N-1 ,, 2ttkn 

X(n) = ]T x(k)w kn with n e [0 ... N-1] and w kn = e N 

k = 0 

-i N-1 

x(n) = — T X(k)w" kn with n e [0 ... N-1] 

N k=0 

In 1965, J.W. Cooley and J.W. Tukey described an 
algorithm allowing to rapidly calculate the Fourier 
5 transform of a series of N complex numbers in an 
article entitled "An .Algorithm for the Machine 
Calculation of Complex Fourier Series", Math. 
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Computation, Vol, 19, 1965, pp. 297-301. This algorithm 
is very interesting when N is a power of two because it 
is then particularly straightforward to implement. This 
algorithm requires p calculation steps where u=log 2 (N) . 
5 It is based on breaking down the series to be 
transformed into two interleaved sub-series. There are 
two kinds of interleaving: time interleaving and 
frequency interleaving. Both kinds of interleaving are 
explained more in detail in the course of the 

10 description with reference to Figs. 1A and IB. 

Specific algorithms have been developed from this 
basic algorithm to deal with the case of real number 
series. The Fourier transformation of a series of 8 
real numbers according to a time interleaving algorithm 

15 and according to a frequency interleaving algorithm is 
illustrated in Figs. 1A and IB. For each time 
interleaving Fourier transform algorithm there is a 
frequency interleaving algorithm, which corresponds to 
a double inversion of the series of transformation 

20 operations, on the one hand, and for each butterfly 
circuit, of the proper transformation operations, on 
the other hand. Whatever the interleaving chosen, the 
transformation method requires three transformation 
steps E0, El, and E2, these steps being implemented 

25 through a set of four butterfly circuits CC, generally 
known as "butterflies" in technical speak. Each 
butterfly, represented in the figure by a point of 
intersection between two columns of numbers, performs 
calculations on two numbers, real or complex. The 

30 symbols R and C respectively identify a real number and 
a complex number. The sequence of real and/or complex 
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numbers after the transformation steps depends on the 
interleaving chosen . 

The time interleaving algorithm is generally 
chosen for calculating the Fourier transform of a 
5 series of real numbers because of the symitietr ical 
distribution of real and complex numbers throughout the 
steps. On the other hand, the frequency interleaving 
algorithm is more suitable for the direct or inverse 
transformation of a series of conjugated complex 
10 numbers. 

When the series to be transformed x(n) is real, 
the Fourier transform verifies the following relation: 
x(n) is real if and only if 

X(n) = X* (-n)=X* (N-n) ; (1) 
15 where * designates the conjugating operation. 

For a series x(n) of N real numbers, the following 
results are inferred from this relation: 

- X(0) and X(N/2) are real; 

- X(n) = X* (N-n) for 1 < n < N/2 -1 

20 Relation (1) highlights the presence of redundant 

information in the X(n) series. 

It should be noted that the transformation method 
is generally implemented by a single set of 
butterflies, the operating mode of which is modified as 

25 the transformation goes along. At each change of 
operating mode, results are stored in a memory having N 
storage locations, the output samples of a butterfly 
replacing the corresponding input samples of the same 
rank in the memory. This method of applying the 

30 algorithm is generally known as the "in place" method. 
This method has a major advantage: if the elements of 



the x(n) series are processed in the first 
transformation step in bit-reversed order of index n, 
the numbers of the X(n) series are output in the last 
transformation step in ascending order of index n and 
vice versa. 

A known transformation method is shown by way of 
example in Fig. 2. This method performs the Fourier 
transformation- of a real x(n) series according to a 
complex time interleaving algorithm. In this example, 
the x(n) series to be transformed comprises sixteen 
real samples, x(0) to x(15). The transformation method 
comprises four transformation steps E p with 0<p<3 . The 
samples of the x(n) series are shown in the first 
transformation step in bit-reversed order of their 
index n. 

At this stage of the explanations, the terms used 
in the course of the description should be defined. The 
rank of a sample is taken from the position it occupies 
in the series of samples to which it belongs. The index 
of a sample then corresponds to the starting rank of 
this sample. 

The intermediate results obtained in the various 
transformation steps are represented by the series 
A(n) , B(n), and C(n). The samples of the series x(n), 
A(n), B(n), C(n), and X(n) are stored in double storage 
locations, one storage location being reserved for the 
real portion of the sample and the other location being 
reserved for the imaginary portion thereof. A R (n) and 
Ai (n) respectively designate the real portion and the 
imaginary portion of the index n sample of the A(n) 
series. Butterflies are represented in the figure by 
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points of intersection between columns of storage 
locations. Each butterfly is assigned a coefficient W s 
symbolized in Fig. 2 by a pair. of coordinates A/B where 
A and B respectively designate the real portion and the 
5 imaginary portion of the coefficient W s . Coordinates 
1/0 and 0/-1 are respectively assigned to coefficients 
W°=l and W N / 4 =W 4 =-j . For the sake of clarity and in 
order to simplify their formulation, the remaining 
coefficients W s have been represented by the following 
10 pairs: 

W 1 -> 2 / -4 W 5 -> -4 / -2 

W 2 -> 3 / -3 W 6 -> -3 / -3 

W 3 -> 4 / -2 W 7 -> -2 / -4 

These pairs of coordinates are graphically 

15 represented in Fig. 3. In fact, coordinates A and B 
respectively represent a cosine value and a sine value. 
This coefficient W s participates in the calculation 
performed by the butterfly. Furthermore, the 
butterflies are distributed at each transformation step 

20 among N/2 P+1 calculation blocks, each calculation block 
comprising 2 P butterflies. In the course of the 
description, the parameter q designates the rank of the 
calculation blocks within the same transformation step; 
q is included between 0 and (N/2 P+1 )-1. 

25 In the first transformation step E 0 , the 

butterflies are distributed among eight calculation 
blocks, each comprising a butterfly performing an 
operation on two complex or real samples. If el and e2 
are to designate the samples applied to the inputs of a 

30 butterfly, the latter outputs samples si and S2 defined 
as follows: 




si = el + W s .e2 and s2 = el -(W s .e2) 
where W s is the coefficient assigned to said butterfly. 

For this first transformation step, coefficient 
W°=l is assigned to the eight butterflies. As the 
5 samples x(n) and the coefficient W° are real, the 
samples A(n) obtained at the end of step E 0 are real. 

For the second transformation step, Ei, the 
butterflies are distributed among four calculation 
blocks each comprising two butterflies. Coefficient 

10 W°=l is assigned to the first one of these butterflies; 
thus, the first butterfly of each calculation block 
provides two real samples. The second butterfly of the 
calculation blocks is associated with coefficient 
W N/4 =W 4 =-j and generates two conjugated complex samples. 

15 The output samples obtained at the end of step Ei are 
designated by the series B (n) . 

For the third transformation step, E 2 , the 
butterflies are distributed among two calculation 
blocks each comprising four butterflies respectively 

20 associated with coefficients W°, W 2 , W 4 , and W 6 . The 
output samples in step E 2 are designated by the series 
C(n). Finally, for the fourth transformation step, E 3 , 
a single calculation block comprising eight butterflies 
respectively associated with coefficients W°, W 1 , W 2 , 

25 W 3 , W 4 , W 5 , W 6 , and W 7 , is provided. This transformation 
step generates the transformed series X (n) . 

Given relation (1), the X(n) series comprises, on 
the one hand, real samples, X(0) and X(8), and on the 
other hand, complex samples, X(l) to X(7) and X(9) to 

30 X(15), samples X(15) to X(9) respectively being the 
conjugates of samples X(l) to X(7). The X (n) series 
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therefore contains redundant information. The storage 
locations outlined in thick stroke in Fig. 2 designate 
the storage locations enclosing the conjugated values 
of the complex samples contained in the storage 
5 locations associated therewith by an arrow. The 
intermediate results series B(n) and C(n) also contain 
redundant information . 

It is then possible to delete this redundant 
information in order to reduce by half the size of the 
10 sample storage memory as well as the number of 
butterflies . 

However, removing redundant information stored in 
the storage locations outline in bold strokes in Fig. 2 
implies the complete reorganization of the 
15 transformation steps of Fig. 2. Thus reorganizing the 
transformation has the effect of modifying the output 
sequence of samples X(n). 

The problem therefore consists in reducing the 
size of the storage memory and the number of 
20 butterflies while maintaining the output sequence of 
the X (n) samples. It is the object of the invention to 
offer a method of calculating the fast Fourier 
transform or the inverse fast Fourier transform of a 
series of N real samples x(n), with N power of 2, 
25 operating according to a time interleaving algorithm, 
which provides the series of samples X(n) in ascending 
order of index n and uses limited calculation and 
storage means . 

For this purpose, the object of the invention is a 
method of calculating the fast Fourier transform or the 
inverse fast Fourier transform of a digital signal 
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defined by a series of N real starting samples x(n), 
with N power of two and n e [0..N-1], comprising 
successive transformation steps for transforming input 
samples into output samples, all the transformation 
steps being performed by means of a single set of 
butterflies with several inputs and several outputs, 
the operating mode of which is modified selectively in 
each transformation step, the input and output samples 
of each transformation step being stored in a storage 
memory, a series of N output samples y(n) 
representative of the fast Fourier transform or the 
inverse fast Fourier transform of the starting samples 
x(n) being provided in the last transformation step, 

characterized in that output samples y(n) are 

real, 

and in that the output samples of a butterfly 
replace the corresponding input samples of the same 
rank in the storage memory, so that, if the starting 
samples x(n) processed in the first transformation step 
are classified in bit-reversed order of their index n, 
output samples y(n) are provided in the last 
transformation step in ascending order of index n, 
these output samples being defined by the following 
relations : 

y(0) = Re[X(0) ] 

y(n) = Re [X ( (n+1) /2) ] for n being odd and 

different from N-l 

y(n) = Im[X(n/2)] for n being even and 

different from 0 

y (N-l) - Re [X (N/2) ] 
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where the X (n) samples, with n e [0..N-1] designate 
the complex samples of the series corresponding to the 
fast or inverse fast Fourier transform of the series of 
starting samples x (n) . 
5 For methods operating according to a frequency 

interleaving algorithm, the invention also relates to a 
method of calculating the fast Fourier transform or the 
inverse fast Fourier transform of a digital signal 
defined by a series of N complex samples X(n) 
10 conjugated by pairs represented by a series of N real 
starting samples y(n), with N power of two and n e 
[0..N-1], the starting samples y(n) being defined as 
follows : 

y(0) = Re[X(0)] 

15 y(n) = Re [X ( (n+1) /2) ] for n being odd and 

different from N-l 
y(n) = Im[X(n/2)] for n being even and 

different from 0 

y (N-l) = Re[X(N/2) ] 

20 this calculation method comprising successive 
transformation steps for transforming input samples 
into output samples, a series of N output samples x (n) 
representative of this fast or inverse fast Fourier 
transform being provided in the last transformation 

25 step, all the transformation steps being performed by 
means of a single set of butterflies with several 
inputs and several outputs, the operating mode of which 
is modified selectively in each transformation step, 
the input and output samples of each transformation 

30 step being stored in a storage memory, 
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characterized in that output samples x(n) are 

real, 

and in that the output samples of a butterfly 
replace the corresponding input samples of the same 
rank in the storage memory, so that, if the starting 
samples y(n) processed in the first transformation step 
are classified in ascending order of index n, the 
output samples x(n) are output in the last 
transformation step in bit-reversed order of index n. 

The inventive calculation methods perform 
operations on real samples and consequently use limited 
calculation and storage means in comparison with the 
method shown in Fig. 2. 

According to another feature of the invention, in 
each transformation step, the butterflies transform 
input sample pairs, the ranks of the input samples of 
the same pair within the series of input samples of 
said transformation step being symmetrical with respect 
to a center between the end rank values of the input 
samples transformed by said butterfly. Input samples 
processed by the same butterfly are thus symmetrically 
linked together by pairs. The result is simplified 
handling of sample addressing. 

According to another aspect of the invention, the 
method preferably comprises p-1 transformation steps E p 
with yt=log 2 (N) and p e [0..p-2]. 



Other features and advantages of the invention 
will be apparent from reading the following detailed 
description, which is made with reference to the 
appended drawings, where: 




Figs. 1A and IB, already described, 
respectively represent a Fourier transformation of 
eight real numbers according to a time interleaving 
algorithm and according to a frequency interleaving 
algorithm; 

Fig. 2, already described, illustrates the 
transformation of a series of 16 real numbers into a 
series of 16 complex numbers according to a complex 
time interleaving algorithm; 

Fig. 3, already described, graphically 
represents the mapping of coefficients W s and 
coordinate pairs A/B; 

Fig. 4 represents modifications applied to part 
of the transformation of Fig. 2; 

Fig. 5 illustrates a modified transformation 
only processing real numbers; 

Fig. 6 is a representation of the method of 
calculating the fast Fourier transform according to the 
invention; 

Figs. 7A and 7B respectively illustrate 
permutations performed on peripheral butterflies and on 
internal butterflies with odd-numbered rank of the 
transformation of Fig. 5; 

Fig. 8 represents an embodiment of the 
transformation method according to the invention, 
comprising y transformation steps; 

Fig. 9 illustrates a grouping of peripheral 
butterflies according to a first embodiment of a 
transformation method comprising y-1 transformation 
steps; 
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Fig. 10 represents a first embodiment of a 
transformation method comprising p-1 transformation 
steps; 

Fig. 11 represents a butterfly design relating 
5 to the transformation method illustrated in Fig. 10; 

- . Fig. 12 represents an alternative of the 
embodiment of Fig. 10; 

Fig. 13 represents a butterfly design relating 
to the embodiment of Fig. 12; 
10 - Fig. 14 illustrates a grouping of peripheral 

butterflies according to a second embodiment of a 
transformation method comprising transformation 
steps; 

Fig. 15 represents a second embodiment of a 
15 transformation method comprising y-1 transformation 
steps, with y being even; 

Fig. 16 represents an alternative of the 
preceding embodiment, with p being odd; 

Fig. 17 represents a butterfly design relating 
20 to the embodiments shown in Figs. 15 and 16; 

Fig. 18 represents a third embodiment of a 
transformation method comprising yi-1 transformation 
steps; 

Fig. 19 represents a butterfly design relating 
25 to the embodiment of Fig. 18; 

Fig. 20 represents the addresses that are 
associated with the various butterflies implemented in 
the embodiment shown in Fig. 12; 

Fig. 21 represents the addresses that are 
30 associated with part of the samples of a transformation 
method processing a series of 32 real samples; 
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- Fig. 22 represents a sample embodiment of a 
transformation method operating according to a 
frequency interleaving algorithm; 

- Fig. 23 represents a butterfly design relating 
5 to the embodiment of Fig. 22; 

Fig. 24 represents the addresses that are 
associated with the various butterflies of the 
embodiment illustrated in Fig. 22. 
G^s According to the invention, only one part of the 

10 X(n) samples is calculated, the other part of the 
samples being redundant. E.g., the calculation could be 
limited to the X (n) numbers for 0<n<N/2. However, a 
more adequate solution is calculating the X(n) numbers 
for 0<n<N/4 and N/2<n<3N/4. The latter solution is 

15 preferable because it only involves the first (N/4)+l 
butterflies for calculating X(n), as the last (N/4)-l 
butterflies can be removed. The storage locations thus 
released can be used for storing the real portion or 
the imaginary portion of the remaining X(n) numbers. 

20 The size of the storage memory can thus be halved if 
the size of the storage locations is limited to storing 
a real number instead of a complex number. The real 
portion of the complex X(n) number is stored in the 
storage location that has been assigned thereto 

25 initially whereas the imaginary portion is stored in 
the storage location initially assigned to the number 
X(N-n) . The same operation can be performed for the 
series of intermediate results B(n) and C(n). The 
transformation method of the invention is limited to 

30 calculating N real samples y(n) instead of N complex 
samples. The real samples y(n) are defined as follows: 
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y (0) = Re[X(0) ] 

y(n) = Re [X ( (n+1) /2) ] for n being odd and 

different from N-l 
y(n) = Im[X(n/2)] for n being even and 

5 different from 0 

y (N-l) = Re [X (N/2) ] . 

Removing redundant information and reorganizing 
storage locations deeply modifies the diagram of 
Fig. 2. Fig. 4 illustrates the reorganizations 

10 performed on the part relating to the calculation of 
samples Ad), AO), A(5) , A(7), B(l), B(3), B{5), B(7), 
C(l) , C(3), C(5), and C(7). Redundant samples to be 
removed are BO), B(7), C(5), and C(7). When the size 
of the storage locations has been reduced and the 

15 samples remaining in theses storage locations have been 
reorganized, the butterflies perform calculations on 
real numbers. The butterflies to which the pair 0/-1 is 
assigned perform calculations on two real numbers. In 
practice, they copy onto their first output the number 

20 present at their first input and multiply by -1 the 
number present at their second input and provide it to 
their second output. The butterflies to which the pair 
1/0 is assigned perform an addition and a subtraction 
on two real numbers. Finally, the other butterflies 

25 perform operations on four real numbers. 

The transformation thus rearranged is illustrated 
in Fig. 5. In this figure, the butterflies associated 
with the pairs 1/0 and 0/-1 corresponding to the 
coefficients W° and W 4 are connected with storage 

30 locations by thick strokes. This figure shows that the 
reorganization of the transformation method steps 
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modifies the output sequence of the y(n) samples and 
therefore the output sequence of the X(n) series. 
Furthermore, this transformation method no longer has a 
specific symmetry allowing to link by pairs the samples 
5 processed by the same butterfly. The result is a very 
complicated address management of the samples to be 
applied to the butterfly inputs. 

Fig. 6 shows the steps of a method of calculating 
the fast Fourier transform or the inverse fast Fourier 

10 transform of a series of N real numbers x (n) , with N 
power of 2, operating according to a time interleaving 
algorithm. It mainly comprises transformation steps 2 
for transforming N starting samples x(n) classified in 
bit-reversed order of their index n into real output 

15 samples y(n) representative of this Fourier transform 
classified in ascending order of index n. 

Advantageously, it comprises a preliminary step 1 
for ranking the N real starting samples x(n) to be 
transformed in the bit-reversed order of their index n 

20 if the samples x (n) are not already in this sequence 
and a final step 3 for generating the N complex samples 
X(n) corresponding to the fast Fourier transform of the 
starting samples x(n) from the N real samples y(n) 
obtained at the end of the transformation steps. 

25 The methods that will be detailed in the course of 

the description will be more in particular for 
calculating the fast Fourier transform of a real 
series. Also, the coefficients W s assigned to the 
butterflies for implementing the inventive method will 

30 be of the type e " jt2 * s/N) with s>0 . For calculating the 
inverse fast Fourier transform the calculation method 



is the same, however, the coefficient is of the type 

e j(2ns/N) with s > 0- 

In order to obtain at the same time y(n) samples 
sorted in ascending order of index n and symmetry of 
calculation, according to the invention, it is 
suggested to modify the calculations performed by the 
butterflies of the odd rank q calculation blocks of the 
transformation illustrated in Fig. 5 in accordance with 
the diagrams of Figs. 7A and 7B. 

As the butterflies associated with coefficient 1/0 
(Fig. 7A) of the odd rank calculation blocks are 
involved, provision is made for permutation both 
outputs of the butterfly and multiplying by -1 the 
result provided at the second butterfly output. 

As the four input butterflies (Fig. 7B) are 
involved, provision is made, for permutation the first 
two outputs with the last two ones. 

This method is applied to the whole transformation 
and then a method is obtained providing at the output 
y(n) samples in ascending order of index n. This method 
is illustrated in Fig. 8. The symbol^ placed above the 
calculation blocks designates the calculation blocks 
wherein butterflies have been modified, i.e. odd rank 
calculation blocks. Due to the symmetry of calculation 
of the transformation method, the intersecting points 
representing the butterflies are superposed inside each 
calculation block. 

The butterflies to which the coefficient 1/0 is 
assigned are called peripheral butterflies because they 
perform calculations on samples arranged at the ends of 
the calculation block. The other butterflies are called 
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internal butterflies. It should be noted that in each 
transformation step, not all the samples are processed 
always and that unprocessed samples are kept in their 
storage locations to be processed in subsequent steps, 
5 or else produced as output if they are already in their 
final shape. 

The transformation method thus modified provides 
y(n) samples in ascending order of index n and, in each 
transformation step, has symmetry of calculation 
10 facilitating the addressing of the samples to be 
processed. 

According to another aspect of the invention, the 
calculation method advantageously comprises 
transformation steps. Several embodiments derived from 
15 that of Fig. 8 and comprising transformation steps 

are therefore shown in the course of the description. 
One butterfly design is associated with each of these 
embodiments. 

All of these embodiments have the following 
20 features in common: 

in each transformation step, provision is made 
for N/2 P+2 calculation blocks and each calculation block 
comprises an peripheral butterfly and/or 2 P -1 internal 
butterflies; all butterflies, be they peripheral or 
25 internal, perform calculations on four real samples; 

the ranks of the samples processed by the same 
butterfly are defined as follows: if, in transformation 
step Ep, a peripheral butterfly belonging to the rank a 
calculation block is taken into consideration, it 
30 transforms the input samples of rank 2 p+2 a, 2 p+2 a+2 p+1 -l , 
2 p+2 a+2 p+1 , 2 p+2 a+2 p+2 -l into output samples of the same 
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rank, and, if an internal rank x butterfly in a rank a 
calculation block is taken into consideration in step 
Ep, it transforms the input samples of rank 2 p+2 a+2x+l, 
2 p+2 a+2x+2, 2 p+2 a+2 p+2 -2x-3, 2 p+2 a+2 p+2 -2x-2 into output 
5 samples of the same rank, with P>1; 

the coefficient assigned to the internal rank x 
butterfly of the rank a calculation block in step Ep is 
equal to W s with 5 = (x+1) . (N/2 P+2 ) . 

In all of these embodiments, the input samples of 

10 each butterfly form sample pairs, the sample ranks of 
the same pair in the series of input samples of a 
transformation step being symmetrical with respect to 
the center value of the end ranks of the input samples 
transformed by said butterfly. This center corresponds 

15 to the value 2 p+2 a+2 p+1 -l/2 . Therefore, one just has to 
know the rank of two of the four samples to be applied 
to the inputs of the butterfly in order to infer 
therefrom the rank of the other two. Addressing these 
samples is thus simplified. This will be explained more 

20 in detail at a later point of the description. 

Thus, according to a first embodiment, the 
neighboring calculation blocks in each transformation 
step are grouped by pairs. The peripheral butterflies 
of the same calculation block are then merged into a 

25 single peripheral butterfly. A sample merging of two 
peripheral butterflies is shown in Fig. 9. This example 
relates to peripheral butterflies associated with 
samples x(0), x(8), x(4), and x(12). 

Furthermore, as the transformation steps do not 

30 process all samples each time, certain calculations can 




be anticipated. E.g., calculating samples C R (1), Ci(l), 
C R (3), and Ci(3) can be done in the second 
transformation step. The result is the diagram of Fig. 
10 showing a first embodiment of the transformation 
5 circuit wherein the method only comprises p-1 
transformation steps. In transformation step E p , each 
calculation block has a peripheral butterfly and 2 P -1 
internal butterflies. It can be considered that this 
method only comprises 3 transformation steps, the 
10 fourth step being limited to performing an addition and 
a subtraction. This addition and this subtraction are 
preferably performed during the final step, and in 
order to limit the number of transformation steps. 

A butterfly design associated with the embodiment 
15 of Fig. 10 is represented in Fig. 11. It comprises: 

four inputs for receiving input samples el,. e2, 
e3, e4, and four outputs for providing output samples 
si, s2, s3, s4, and 

three additional, respectively primary mode MP, 
20 permutation PERM, and coefficient COEF, inputs. 

This butterfly is responsible for selectively 
applying to input samples el, e2, e3, and e4, various 
transformation operations each determined by the values 
assigned to primary mode, permutation signals and to 
25 coefficient W s admitted at the corresponding additional 
inputs . 

The primary mode signal is 0 for a peripheral 
butterfly and 1 for an internal butterfly. When the 
permutation signal is a 1, the output samples si and s2 
30 of the butterfly are swapped with output samples s3 and 
s4. This permutation is only possible if the butterfly 
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is an internal one. Finally, the coefficient W s 
associated with the butterfly is applied to the 
coefficient input COEF. 

Thus, if the complex coefficient W s =A+j.B is 
5 applied to the coefficient entry of the butterfly, the 
latter provides the following output samples si, s2, 
s3, and s4 

1) if the primary mode signal is 0: 
si = el + e2 

10 s2 = el - e2 

s3 = e4 - e3 
s4 = e3 + e4 

2) if the primary mode signal is 1 and the 
permutation signal is 0: 

15 si = el + A.e3 - B.e4 

s2 = e2 + B.e3 + A.e4 
s3 = el - A.e3 + B.e4 
s4 = -e2 + B.e3 + A.e4 

3) if the primary signal is 1 and the permutation 
20 signal is 1 : 

si = el - A.e3 + B.e4 

s2 = -e2 + B.e3 + A.e4 

s3 = el + A.e3 -B.e4 

s4 = e2 + B.e3 + A.e4 
25 According to an alternative embodiment, provision 

can be made for addition and subtraction to be 
performed inside the peripheral butterfly of the last 
transformation step. This embodiment is shown in Fig. 
12. For this purpose, the corresponding butterfly 
30 design has a fourth additional input called secondary 
mode input MS to which a secondary mode signal is 
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applied. This signal is 1 for the peripheral butterfly 
of the last transformation step, otherwise it is 0. 
This design is illustrated in Fig. 13. This design has 
an additional operating mode in comparison with the 
5 preceding one; thus, when the primary mode signal is 0 
and the secondary mode signal is 1, the output obtained 
is : 

si = el + e2 + e3 + e4 ' 
s2 = el - e2 
10 s3 = e4 - e3 

s4 = (el + e2) - (e3 + e4) 

According to a second embodiment derived from the 
diagram of Fig. 8, the neighboring peripheral 
butterflies in even index p transformation steps E p are 

15 grouped by pairs and are merged with the peripheral 
butterfly of the second odd index step in order to form 
a new peripheral butterfly at the odd index step. This 
grouping is illustrated in Fig. 14 in an example. In 
this example, the peripheral butterflies of the first 

20 transformation step processing samples x(0), x(8), 
x(4), and x(12) are merged with the peripheral 
butterfly of the second step processing samples A R (0) 
and A R (2). Also, the peripheral butterflies of the 
first transformation step processing samples x(2), 

25 x(10), x(6), and x(14) are merged with the peripheral 
butterfly of the second step processing samples A R (4) 
and A R (6). The two butterflies obtained are different 
in that the second one performs in addition a 
permutation between the first and second outputs. If 

30 this grouping is applied to the whole transformation 
illustrated in Fig. 8, the result is that step E 0 no 
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longer uses any butterflies and can be removed. The 
resulting transformation method is illustrated in Fig. 
15. 

However, two cases should be distinguished for 
5 this transformation method: the case where N is an even 
power of two (p. being even) and the case where N is an 
odd power of two (p being odd) . 

In case ]i is even, there is an even number of 
transformation steps in the embodiment shown in Fig. 8 
10 and grouping the peripheral butterflies of even index 
steps with those of the following odd index steps is no 
problem. This case corresponds to the diagram in Fig. 
15. 

In case p is odd, the peripheral butterfly of the 

15 last even index step cannot be grouped with other 
peripheral butterflies. Therefore, a specific operating 
mode should be provided for this case. This case is 
illustrated in Fig. 16, this figure representing the 
transformation of a series of eight real samples (u=3) . 

20 The peripheral butterfly of the last transformation 
step of this circuit could not be merged with other 
peripheral butterflies . 

The butterfly design associated with this second 
embodiment is illustrated in Fig. 17; it differs from 

25 the preceding design in that the secondary mode signal 
is 1 when a peripheral butterfly for implementing the 
last step is involved and ]i is odd, and in that 
permutation applies to all the butterflies of the even 
rank calculation blocks. 

30 The calculations performed by the butterfly are 

also different and are defined as follows: 
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1) if primary mode, secondary mode and permutation 
signals are 0 : 

si = el + e2 + e3 + e4 
s2 = el - e2 
5 s3 = e4 - e3 

s4 = (el + e2) - (e3 + e4) 

2) if the primary mode signal is 0 and the 
secondary mode signal is 1: 

si = el + e4 
10 s2 = e2 

s3 = e3 
s4 = el - e4 

3) if the primary mode signal is 0 and the 
permutation signal is 1: 

15 si = (e3 + e4) - (el + e2) 

s2 = el - e2 
s3 = e4 - e3 
s4 = el + e2 + e3 + e4 

4) if the primary mode signal is 1 and the 
20 permutation signal is 0: 

si = el + A.e3 - B.e4 

s2 = e2 + B.e3 + A.e4 

s3 = el - A.e3 + B.e4 

s4 = -e2 + B.e3 + A.e4 
25 5) if the primary mode signal is 1 and the 

permutation signal is 1: 

si = el - A.e3 + B.e4 

s2 = -e2 + B.e3 + A.e4 

s3 = el + A.e3 - B.e4 
30 s4 = e2 + B.e3 + A.e4 



In case y is odd, it is also possible on the one 
hand to provide for the peripheral butterflies to be 
grouped for implementing the first transformation step 
in the same way as in the first embodiment/ and on the 
other hand, for the butterflies of the other steps to 
be grouped as in the third embodiment. Grouping 
butterflies from the second step on is then performed 
by assigning an even index to the first transformation 
step- These groupings are represented in Fig. 18. 

The butterfly design corresponding to this 
embodiment is represented in Fig. 19. The secondary 
mode signal is 1 for a peripheral butterfly 
implementing the first transformation step of the 
circuit and if y is even. The calculations performed by 
this butterfly are the following ones: 

1) if the primary mode signal is 0 and the 
secondary mode signal is 1: 

si = el + e2 
s2 = el - e2 
s3 = e4 - e3 
s4 = e3 + e4 . 

2) if primary mode, secondary mode and permutation 
signals are 0: 

si = el + e2 + e3 + e4 
s2 = el - e2 
s3 = e4 - e3 

s4 = (el + e2) - (e3 + e4) 

3) if the primary mode and secondary mode signals 
are 0 and the permutation signal is 1: 

si = (e3 + e4) - (el + e2) 
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s3 = e4 - e3 
s4 = el + e2 + e3 + e4 

4) if the primary mode 
permutation signal is 0: 

5 si = el + A.e3 - B.e4 

s2 = e2 + B.e3 + A.e4 
s3 = el - A.e3 + B.e4 
s4 = -e2 + B.e3 + A.e4 

5) if the primary mode 
10 permutation signal is 1: 

si = el - A.e3 + B.e4 
s2 = -e2 + B.e3 + A.e4 
s3 = el + A.e3 - B.e4 
s4 = e2 + B.e3 + A.e4 
15 In all the embodiments described before, the ranks 

of samples processed by the same butterfly are 
symmetrical by pairs with respect to a center value. 
One just has to know the rank of the first two input 
samples of the butterfly to infer the other two ones 
20 therefrom by symmetry. If the input and output samples 
are saved in storage locations the address of which 
corresponds to the rank of these samples, addressing 
the latter will be simplified. 

Indeed, all that needs to be done then is to 
25 generate two addresses per butterfly, as the other two 
ones can be inferred by symmetry. Furthermore,- it 
should be noted that the addresses of the input samples 
and those of the corresponding output samples are the 
same because the transformation is applied according to 
30 an in place method. 
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The addresses associated with the various 
butterflies relating to the transformation method of 
Fig. 12 are shown in Fig. 20. The address of a sample 
is taken to be equal to the rank of this sample in the 
5 sample series to which it belongs. The series to be 
transformed in the example of Fig. 12 comprises 16 
samples. Therefore, 16 addresses have to be produced, 
with four-bit addressing. For a series of N samples, 
log2(N)-bit addressing is performed. 

10 Each transformation step of the method of Fig. 12 

is implemented by four butterflies each processing four 
real input samples. The binary addresses of the samples 
to be processed in each step are therefore distributed 
among four groups of four addresses. The address groups 

15 relating to one peripheral butterfly are contained in 
bold line boxes and the address groups relating to an 
internal butterfly are contained in thin line boxes. 
Furthermore, the address groups relating to the same 
calculation block are grouped in dotted line boxes. 

20 The first peripheral butterfly in step E 0 

processes the first four samples of the x(n) series 
resulting from the preliminary classification step. The 
corresponding addresses to be generated for this 
butterfly are therefore 0000, 0001, 0010, and 0011. 

25 Also, considering the first peripheral butterfly of 
step Ei, it processes the samples of rank 0, 3, 4, and 
7 of the A(n) series. The corresponding addresses to be 
generated for this butterfly are therefore 0000, 0011, 
0100, and 0111. 

30 The (u-2-p) high-order bits of the addresses 

indicate rank q of the calculation block to which the 
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address is related. Thus, the two high-order bits of 
the addresses generated for the first calculation block 
of step E 0 are 00. It should also be noted that, as the 
last step only uses a single calculation block, the 
5 addresses generated for this calculation block have no 
bit relating to the rank of this block (y-p-2=0) . 

For each butterfly, provision is made for 
generating only two binary addresses; the other two 
ones are obtained by inverting the (p+2) low-order bits 

10 of the generated addresses- E.g., considering the 
addresses of the first peripheral butterfly in step Ei, 
only addresses 0000 and 0011 are generated and 
addresses 0111 and 0100 are obtained by inverting the 3 
low-order bits of the generated addresses. 

15 A first and a second address per butterfly are 

then produced, these addresses being consecutive for an 
internal butterfly. When a peripheral butterfly is 
involved, the p+2 low-order bits of the first address 
are equal to 0, and the p+2 low-order bits of the 

20 second address form a number equal to 2 P+1 -1. 

Regarding the coefficients W s to be applied to the 
coefficient inputs COEF of the internal butterflies, 
they are stored in a memory of the calculation circuit. 
Only N/4 coefficient values are required for 

25 calculating the Fourier transform, yi-2-bit addressing 
is performed for these coefficients. In the embodiments 
described previously where N=16, only the pairs 1/0, 
2/-4, 3/-3, and 4/-2 are used which correspond to 
coefficients W°, W 1 , W 2 , W 3 . 

30 According to the invention, the address of these 

four coefficients must therefore be known. The address 
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associated with each coefficient W s is chosen to be 
equal to the value of the power s. Consequently, the 
addresses of coefficients W°, W 1 , W 2 , W 3 are 00, 01, 10, 
11, respectively. 
5 So as not to be obliged to generate these 

addresses, according to the invention, the addresses 
produced for addressing the samples are used. The 
address of the coefficient which is assigned to a 
butterfly is included in the second address produced 
10 corresponding to the greater one of both addresses. 

However, we must distinguish between three cases: 

a) when p+l=p-2, the coefficient address 
corresponds to the number formed by the p+1 low-order 
bits of the second address generated for this internal 

15 butterfly. This is the case for the second step (p=l) 
in the example of Fig. 20. The 2 low-order bits of the 
second address are 10 and therefore designate 
coefficient W 2 . 

b) if p+l>yi-2, the coefficient address corresponds 
20 to the number formed by the p+1 low-order bits of the 

second address generated for this internal butterfly, 
minus its ]i-p-l low-order bits. This is the case for 
the third step (p=2) in the example of Fig. 20. The 3 
low-order bits of the second address generated for the 

25 first internal butterfly are 010. When the last bit 

p-l=l) of this number is taken away, the number 01 is 
obtained relating to coefficient W 1 . This case always 
corresponds to the last step of a transformation method 
comprising p-1 transformation steps. 

30 c) if p+Kp-2, the coefficient address corresponds 

to the number formed by the p+1 low-order bits of the 




second address generated for this internal butterfly, 
followed by y-p-3 zero bits at the end of the number. 
This case is illustrated in Fig. 21. This figure 
represents the addresses relating to an internal 
5 butterfly for implementing the second transformation 
step (p=D of a transformation method designed for 
processing a series of 32 real samples (N=32 and y=5) . 
This butterfly transforms the samples of rank 1, 2, 5, 
and 6 of the series of samples obtained at the end of 

10 the first step. The two low-order bits of the second 
generated address are 10 and if a zero (y-p-3=l) is 
added at the end of the number, the number 100 is 
obtained designating coefficient W 4 . 

Thus, both addresses produced by the address 

15 generator for one butterfly is used to address both 
four samples to be processed and the coefficient 
relating to the butterfly. 

Preferably, even address and odd address samples 
will be stored in two separate memories. Thus, it will 

20 be possible to read two input samples simultaneously, 
and it will be possible to write the resulting output 
samples simultaneously, which means saving process time 
for the series to be transformed. 

As mentioned before, for each calculation method 

25 operating according to a time interleaving algorithm 
there is a corresponding method operating according to 
a frequency interleaving algorithm. All that has to be 
done to obtain it is, on the one hand, to invert the 
series of transformation operations of the 

30 corresponding time interleaving method, and on the 
other hand, for each butterfly, to invert the 
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transformation operations themselves with respect to 
those of the corresponding method. 

In addition, the invention also relates to a 
method of calculating the fast Fourier transform or the 
5 inverse fast Fourier transform of a series of N complex 
samples X (n) conjugated by pairs, with N power of 2, 
operating according to a frequency interleaving 
algorithm. The series of N complex samples X(n) is 
represented by a series of N real samples y(n) defined 
10 as follows: 

y (0) = Re[X(0) ] 

y(n) = Re [X ( (n+1) /2) ] for n being odd and 

different from N-l 
y(n) = Im[X(n/2)] for n being even and 

15 ^ different from 0 

y (N-l) = Re[X (N/2) ] 

According to the invention, this method 
substantially comprises transformation steps for 
transforming input samples into output samples. Real 

20 samples y(n) are processed in a first transformation 
step and the last step provides a series of N output 
samples x(n) representative of the fast or the inverse 
fast Fourier transform of the sample series X(n). As 
for the time interleaving circuits, each transformation 

25 step is implemented by a set of butterflies with 
several inputs and several outputs. As the 
transformation is done according to an in place method, 
all the steps are performed by means of a single set of 
butterflies, the operating mode of which is modified in 

30 each transformation step. In each transformation step, 
input and output samples are stored in a storage 



o gi 5 b ± e 7' e . i. o ± y o is: 



31 



memory. When they have been transformed, the output 
samples of the same butterfly replace the corresponding 
input samples of the same rank in the storage memory. 

According to the invention, if the samples y(n) 
5 introduced in the first transformation step are 
classified in ascending order of index n, the output 
samples x (n) are provided in the last transformation 
step in bit-reversed order of index n. The output 
sequence of the x(n) samples can then be modified by a 
10 final step so as to classify them in ascending order of 
index n. 

An embodiment of such a transformation is 
represented in Fig. 22. It is inferred from the 
embodiment of Fig. 12 by inverting the functional 

15 arrangement of the transformation steps of Fig. 12 
(mirror image of what it is for time interleaving) . 
This embodiment allows the calculation of the inverse 
fast Fourier transform of a series of 16 real samples 
y(n) representative of a series of 16 complex samples 

20 X(n) conjugated by pairs. Coefficients W s are therefore 
of the type e j(27ls/N) . 

This embodiment comprises three transformation 
steps E p with 0<p<2 . In each transformation step E p , 
the butterflies are henceforth distributed among 2 P 

25 calculation blocks, these calculation blocks being 
sorted in each step according to an ascending rank q 
from 0 to 2 P -1. Each calculation block has a peripheral 
butterfly and N/2 p+2 -l internal butterflies. The y(n) 
samples are applied in the sequence of index n in the 

30 first step. 
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The ranks of the samples processed by the same 
butterfly are defined as follows: in step Ep, 
considering a peripheral butterfly belonging to the 
rank a calculation block, it transforms the input 
samples of rank 2^ p ct, 2^ p a+2^ p " 1 -l , 2^ p a+2^ p -\ 2^ p a+2^ 
p -l into output samples of the same rank, and, in step 
Ep, considering an internal rank x butterfly in the 
rank a calculation block, it transforms the input 
samples of rank 2 p - p a+2x+l, 2 p - p a+2x+2, 2 M " p a+2 p " p -2T-3 , 2 P " 
p a+2 p ~ p -2x-2 into output samples of the same rank. 
Finally, the coefficient assigned to the rank x 
internal butterfly of the rank a calculation block in 
step Ep is equal to W 5 with 8 = (x+1) ,2 P . 

The real samples x (n) obtained at the end of the 
method are provided in bit-reversed order of index n. 

At the butterflies, the coefficients W s are of the 
type e j(2 * s/N) and inputs and outputs have been inverted 
with respect to the embodiment of Fig. 12. 
^ Consequently, the operations performed by the 
butterflies of this embodiment are different from those 
performed by the butterflies of Fig. 12. A butterfly 
design associated with the embodiment of Fig. 22 is 
represented in Fig. 23. 

Just like the design associated with Fig. 12, it 
has four data inputs and four data outputs as well as 
four additional inputs, respectively primary mode MP, 
secondary mode MS, permutation PERM, and coefficient 
COEF inputs. The primary mode signal is 0 for an 
peripheral butterfly and 1 for an internal butterfly. 
The permutation signal is 0 for even values of rank q 
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and 1 for odd values. Finally, the secondary mode 
signal is 1 if the peripheral butterfly is used for 
implementing the first step, and otherwise 0. 

The calculations performed by this butterfly are 
the following ones (W s =A+j.B): 

1) if the primary mode and secondary mode signals 
are 0 : 

si = (el + e2)/2 
s2 = (el - e2) /2 
s3 = (e4 - e3) /2 
s4 = (e3 + e4) /2 

2) if the primary mode signal is 0 and the 
secondary mode signal is 1 : 

si = [ (el+e4) /2+e2] /2 
s2 = [ (el+e4) /2-e2] /2 
s3 = - [e3- (el-e4) /2] /2 
s4 = [e3+ (el-e4) /2] /2 

3) if the primary mode signal is 1 and the 
permutation signal is 0: 

si = (el+e3) /2 
s2 = (e2-e4) /2 

s3 = [(el-e3).A - (e2+e4).B]/2 
s4 = [-(el-e3).B + (e2+e4).A]/2 

4) if the primary mode signal is 1 and the 
permutation signal is 1: 

si = [(el-e3).A - (e2+e4).B]/2 
s2 = [-(el-e3).B + (e2+e4).A]/2 
s3 = (el+e3)/2 
s4 = (e2-e4)/2 
This design is inferred from the design of Fig. 13 
by inverting values el, e2, e3, e4 and values si, s2, 
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s3, s4, and by replacing B with -B because coefficient 
W s is now of the type e ^ (27Cs/N) . as A and B represent the 
cosine and sine of the same number, they have become 
A 2 +B 2 =l. The expressions of si, s2, s3, s4 are thus 
5 simplified* 

According to the invention, addresses are 
furthermore generated for sample addressing. Provision 
is made for generating two binary addresses of p bits 
per butterfly, each binary address corresponding to the 

10 rank of a butterfly input sample. The addresses of the 
other two samples to be applied to the butterfly inputs 
are obtained by inverting the ]i-p low-order bits of the 
first two addresses. 

In the same way as for the transformation methods 

15 operating according to time interleaving, the two 
binary addresses produced are consecutive for an 
internal butterfly. For a peripheral butterfly, the p-p 
low-order bits of the first generated address are equal 
to 0, and the p-p low-order bits of the second address 

20 form a number equal to N/2 P+1 -1. By way of example, the 
addresses produced for the transformation circuit of 
Fig. 22 are gathered in Fig. 24. 

Advantageously, provision can be made for storing 
even address samples and odd address samples in two 

25 separate memories in order to reduce processing time of 
the transformation operation. 

Finally, the addresses generated for addressing 
the samples are also used for addressing coefficients 
W s . The value of the parameter s is used for addressing 

30 the corresponding coefficient W s . In this embodiment, 
the parameter s is equal to: 



if y-p-l=u-2 , the number formed by the y-p-1 
low-order bits of the second address produced for said 
internal butterfly, 

if jj-p-Kjj-2, the number formed by the p-p-1 
low-order bits of the second address produced for said 
internal butterfly, followed by p-1 zero bits at the 
end of the number, 

if p-p-l>vi-2, the number formed by the p-p-1 
low-order bits of the second address produced for said 
internal butterfly, minus its p+1 low-order bits. This 
case corresponds to the first step (p=0) of the 
transformation methods operating according to a 
frequency interleaving algorithm. 



